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Abstract 

Decoherence in quantum systems which are classically chaotic is studied. The 
Arnold cat map and the quantum kicked rotor are chosen as examples of linear and 
nonlinear chaotic systems. The Feynman- Vernon influence functional formalism is used 
to study the effect of the environment on the system. It is well-known that quantum 
coherence can obliterate many chaotic behavior in the corresponding classical system. 
But interaction with an environment can under general circumstances quickly diminish 
quantum coherence and reenact many classical chaotic behavior. How effective deco- 
herence works to sustain chaos, and how the resultant behavior qualitatively differs 
from the quantum picture depend on the coupling of the system with the environment 
and the spectral density and temperature of the environment. We show how recurrence 
in the quantum cat map is lost and classical ergodicity is recovered due to the effect 
of the environment. Quantum coherence and diffusion suppression are instrumental to 
dynamical localization for the kicked rotor. We show how environment-induced effects 
can destroy this localization. Such effects can also be understood as resulting from ex- 
ternal noises driving the system. Peculiar to decohering chaotic systems is the apparent 
transition from reversible to irreversible dynamics. We show such transitions in the 
quantum cat map and the kicked rotor and distinguish it from apparent irreversibility 
originating from dynamical instability and imprecise measurements. By performing 
a time reversal on and following the quantum kicked rotor dynamics numerically, we 
show how the otherwise reversible quantum dynamics acquires an arrow of time upon 
the introduction of noise or interaction with an environment. 
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1 Introduction and Summary 



1.1 Quantum versus Classical Chaos 

The problem of quantum chaos has been intensively studied in the recent decade |], || 
H |J ||, ||, [7|, ||, H 0, [ll[]. Although the precise criteria for quantum chaos is still not 
well established at this stage, the salient features of a quantized classically- chaotic system 
are better understood than before. In classical dynamics, chaos appears as the result of 
instability caused by nonlinearity or the compactness of the phase space, as manifested in 
the quantum kicked rotor and the Arnold cat map ||, two examples we will discuss in this 
paper. The degree of instability can be measured by the exponential rate of separation 
of initially infinitesimally-close intervals, namely, the Lyapunov exponent. When this local 
instability occurs for the entire phase space, global chaos sets in. [] 

To look for similar phenomena in quantum systems one encounters basic difficulties. 
To begin with, the very concept of trajectories which is used to define classical chaos is 
meaningless in quantum mechanics. The equations of motion in quantum mechanics is linear. 
Seeking nonlinear effects in these linear equations, as well as using concepts of determinacy 
in a theory based on probabilistic interpretations are intrinsically prohibitive. The words 
'quantum chaos' generally refer to possible traces or shadows of chaos in the quantum system 
obtained from quantizing the corresponding classical system which are known to possess 
chaotic behavior. The study of quantum chaos is devoted to finding how the classical notion 
of instability changes when the system is quantized, and how such changes can be expressed in 
the language of quantum mechanics. For example, fingerprints of classical chaos may appear 
as scars in the wavefunction, as fluctuations in the spectrum, or as diffusion localization, etc 

How are these classical and quantal characteristics related to each other in the corre- 
spondence between quantum and classical chaos? There are many criteria of classicality, 
an issue whose recent resurgence of interest is stimulated by developments in many areas of 
physics (see, e.g., fll^l ). Using the uncertainty principle as one criterion, we see immediately 
that there is a fundamental discrepancy between the definition of chaos and quantum uncer- 



tainty JT3[. For systems with conservative dynamics, the initially infinitesimally-separated 
trajectories in phase space will exponentially diverge in some direction and converge in some 
other direction. This will soon become incompatible with the quantum uncertainty principle 
which prevents one to specify details between points in the phase space separated closer than 
the Planck constant h. Therefore, the fact that many classically-chaotic systems produce 
infinitely- folded, Cantor-set stuctures which can continue to arbitrarily small scale due to 
nonlinearity is in conflict with quantum mechanics. 

In classical mechanics, nonlinearity makes the dynamics sensitive to finer scales, leading 
to various fractal structures. However, in quantum mechanics, owing to the linear nature 



^rgodicity thus generated is one of the basic criteria for the validity of equilibrium statistical mechanics. 
Infinite repetition of streching and folding in the phase space may be the cause for the generation of self- 
organized structures in the microscopic world. 
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of Schrodinger's equation, one expects to see limitations to such fine structures. Quantum 
effects are known to smooth out the many- folded trajectories caused by nonlinearity. In this 



respect, quantum effect is similar to the effect of noise on classically chaotic systems [14 



In fact, one can study the scaling property from the classical to the quantum regimes in a 



system when Ti — > as if the system is subject to some external noise |TH] 



1.2 Classicality as an Emergent Behavior of Quantum Open Sys- 
tems 

The above description of quantum-classical correspondence takes the point of view that 
quantum effect is a correction to the underlying classical dynamics, which is the attitude 
taken by many work on this subject using semi-classical approximations. However, this 
is opposite to how nature works: most of us would agree that quantum mechanics is the 
fundamental theory which describes nature, and classical mechanics is only an approximation 
to it. 

How classical dynamics arises from the fundamental principles of quantum mechanics and 
how our ordinary classical experience can be reconciled with the quantum depiction have 
been the basic questions asked in the foundation of quantum mechanics and in quantum 
measurement theory. Although many different explanations exist, the environment-induced 
decoherence point of view seems to be one of simple and practical importance. |T6| . In this 
point of view, the quantum to classical transition is induced by the interaction of a quantum 
system with an environment. The averaged effect from coarse-graining the large number of 
degrees of freedom is the diminuation of quantum coherence and the appearance of diffusion 
and dissipation in the effective dynamics of the system. The decoherence time is defined as 

W=-(^) 2 (1.1) 
7 ox 



where tdis = 7 _1 is the dissipation time scale (7 is the damping constant), X edB = h/^/2ixmkT 
is the thermal de Broglie wavelength, and Sx is the characteristic size of the system (here 
we assume a coordinate coupling x). The decoherence time is usually very short for a bath 
at high temperatures. We refer the readers to recent reviews on this topic |TB . 



This approach has been applied to problems involving quantum decoherence in quan- 
tum measurement theory, mesoscopic physic, quantum cosmology and semiclassical gravity 



Hq| . However, The quantum and classical correspondence of chaotic systems in terms of 
environment-induced decoherence has so far been studied only by a limited number of au- 
thors 0, |T8|, |I9 . 



1.3 Decoherence, Localization and Irreversibility 
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1.3.1 Noise and Localization 



There are many detailed studies of the classical and quantum kicked rotor model ||. A 
particularly interesting feature of quantum nonlinear chaotic systems is the localization of 
wave functions in momentum space due to quantum coherence p2| |. This momentum- 
space localization of the wave function is often compared with Anderson localization [20| of 
electrons in a random potential. Localization in the kicked rotor is considered to occur by 
a similar mechanism p3| . If one views classicality as an emergent behavior of a decohered 
quantum system, then it is of interest to study the effect of an environment on localization. 
Dittrich and Graham [24| studied the kicked rotor in a harmonic oscillator bath, and derived 
a master equation for the open system. Their argument is mainly focused on the effect 
of dissipation induced by the environment. They used a low temperature approximation 
and in the zero temperature limit they claim that the dynamics becomes Markovian. This 
rather unusual behavior is due to the special non-Ohmic environment they used. In general, 
the Markovian regime corresponds only for an ohmic bath at high temperature. Cohen 



and Fishman |25| used the influence functional method |^6| to study the effect of noise 
associated with an Ohmic bath on localization for the QKR and a similar model. They 
calculated explicitly the diffusion constant and the relevant time scales in terms of the noise 
correlation and the nonlinear parameter. On a related problem, Ott, Antonsen and Hanson 
|2T|| first showed numerically that external noise breaks the localization of a wave packet in 
the QKR. Cohen also studied the effect of noise correlations [25] . Naively one does not expect 
correlations to play an essential role for chaotic systems because the memory in such systems 
is lost quickly. However, in the quantal case, long range correlations may alter the situation 
in a complicated way. In fact it is known that the appearance of noise autocorrelation 
depends on the system-environment coupling. 



1.3.2 Irreversibility 

Using a simple linear continuous model, the inverted harmonic oscillator potential, Zurek 
and Paz [^] recently observed that in the presence of noise, the dynamics can change from 
a reversible classical one to an irreversible one. We show that a similar behavior exists in 
the kicked rotor model. For a conserved Hamiltonian chaotic system, volume conservation 
causes one direction in phase space to contract exponentially. Without interaction with an 
environment, the other source of irreversibility intrinsic to chaotic systems arises from the 
limitation of actual measurements. For example, the classical kicked rotor behaves essentially 
irreversibly due to the instability of trajectories. However, in the quantal case, the system 
characterized by the quantum state becomes highly stable in spite of the nonlinearity of the 
Hamiltonian [[31]]. When this system interacts with an environment, there exists a sharp 
transition from a quantum reversible conservative stage to a classical irreversible stage. 

2 The word "localization" here refers to Anderson localization ^0|, not to the establishment of delta- 
functional correlation between, say, the momentum and the coordinate in the realization of the quasi-classical 
state which is sometimes used in the context of decoherence and quantum to classical transition plf . 
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Irreversibility ascribed to by a limitation of measurement will be replaced by irreversibility 
arising from coarse graining the environment. 



1.4 Time Scales of Competing Processes 

One way to gauge the relative importance of the pertinent processes which can influence the 
dynamics of a quantum chaotic system is to compare their characteristic time scales. Let us 
start with cases with no interaction with an environment. There are essentially two different 
time scales involved. One is the Ehrenfest time £g and the other is the relaxation time £r. 
The Ehrenfest time is defined as the time within which the Ehrenfest theorem holds. 

where Sp(0) is the relevant initial (angular) momentum scale. Violation of the Ehrenfest 
theorem in the quantal case arises from the nonlinear terms in the potential which can be 
seen in the evolution (Kramer-Moyal) equation of the Wigner function. 

dW(X,p) 2 . h d d d x n ^ 

- --Hsm(-( — —-——))W{X,p) (1.3) 



dt h y 2 y dpdX dXdp' 

2 l d 2n+l Hd 2n+l W 

h' (2n + l)\dX 2n+1 dp 2n+1 



= {H,w} + ^ =1 (-in-r fn ' 2l) ;; - n.n 

where { } is the Poisson bracket. 

The appearance of localization arises from the discrete spectrum of the Hamiltonian. 
Thus the time it takes for the wave packet to localize is determined by how long it takes for 
the system to recognize the discreteness of the spectrum. Simple argument is given in [p2| : 
Since I represents the effective number of modes scattered in the period [0, 27r], the typical 
spacing is given by Aoj ~ 2n/l. Thus after £r ~ 1/Au ~ I, the system localizes. (See also 
discussion in Sec. 3.3.) 

Upon interaction with a bath, a system effectively decoheres at the decoherence time 
scale tdec- Another time scale to arising from the coarse graining appears which also con- 
tributes to the violation of the Ehrenfest theorem. As discussed in [[Tj|, it determines the 
transition regime from the reversible classical Liouville dynamics to the irreversible dynamics 
as embodied in the Second Law of Thermodynamics. 

From our study, this picture also holds for the kicked rotor. In this case, we see the 
transition from the initial constant-entropy regime to the entropy- increasing regime. Because 
of the compactness of the space, we see entropy does not increase forever but will eventually 
saturate. After tc, the evolution is not unitary. Note that even if the evolution of the Wigner 
function is the same as that of the classical Liouville distribution function before one 
should not regard the system as in a classical state. 



3 Note that, as shown in the examples of 32 , the Ehrenfest theorem is neither necessary nor sufficient to 
define classicality. There are systems which evolve strictly quantum mechanically but the expectation values 
of the canonical variables obey classical equations; and there are models which do not satisfy the theorem 
but their evolution is essentially classical. 
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Dynamical localization is completed at the relaxation time £r. At decoherence time td ec , 
coherence is destroyed up to the localization length. If t dec >> t R , suppression of momentum 
diffusion due to quantum effects always exists and we will never see the classical state. 

It is known that the evolution of the kicked rotor is not time reversible while its quantized 
version is time reversible. This type of quantum stability is also considered to be one of 
the characteristics of quantum chaos. In a quantum system, irreversibility arising from 
limited precision in a measurement now no longer causes serious loss of information. Instead, 
interaction with a bath introduces the irreversibility due to the coarse graining for the quantal 
case. 

In this paper, we study the quantum dynamics of two simple models which possess 
classical chaotic behavior, the Arnold cat map and the kicked rotor. By introducing linear 
coupling with a harmonic oscillator bath assumed to be Ohmic and at high temperature, we 
show how the effective dynamics of a quantum open system reveals the well-known classical 
chaotic behavior. In Section 2, we examine the quantum cat map (QCM) of a system coupled 
with a harmonic bath. The system is known to be chaotic when the corresponding matrix 
for the map is hyperbolic. We use the influential functional method to study the effect of 
the environment on this system. By measuring the linearized entropy we show that the 
decoherence mechanism works more efficiently than the regular case. Namely, the rate of 
decoherence is faster in the chaotic system. Decoherence rate in chaotic systems was also 
studied by Tameshit and Sipe |l7fl. Peculiar to the quantum case is the recurrence behavior 
of physical quantities, resulting from the finiteness of the phase space points in the quantum 
map due to the quantization ( because the phase space is periodic in both the coordinate 
and momentum). We show that interaction with the environment erases the recurrence 
in the hyperbolic map but not in the elliptic map. Thus both cases behave close to the 
correponding classical limit. 

In Section 3, we examine the quantum kicked rotor (QKR) as a prototype of nonlinear 
chaotic systems. Without interaction with a bath, the wave function shows localization 
arising from quantum coherence effects. Loss of coherence due to interactions with an envi- 
ronment shown by the decay of the off-diagonal components of a reduced density matrix is 
responsible for the breaking of localization. The decay rate increases as the noise strength 
associated with the environment and the nonlinear parameter get larger. In Section 4, we ex- 
amine the transition from reversible to irreversible dynamics intrinsic to an unstable system 
due to the interaction with an environment. We show that the same mechanism holds for 
both the cat map and kicked rotor. For both cases, the entropy shows saturation possibly due 
to the bounded nature of the phase space. Furthermore, we perform time-reversal numeri- 
cally and show how the interaction with an environment changes the nature of irreversability. 
Details of results can be found at the end of each section. 



2 Decoherence in a Linear Map 
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2.1 Quantum Cat Map 

The cat map is a linear area-preserving map T on a two-torus in phase space by identifying the 
boundaries of the interval [0, 27r] in both the coordinate Q and the momentum P directions 
From time step j to j + 1, it is given by 



(&)-(: s)(2M2) 

where detT = 1 guarantees area preservation. The degree of chaos depends on the choice 
of T. The eigenvalues of T are either both real or both imaginary. In the latter case, T is 
elliptic, the motion becomes periodic and no sensitive dependence on the initial condition is 
observed. When T is hyperbolic, the motion becomes chaotic. 

Quantized cat map is studied in detail by Hannay and Berry f54]| . Due to the periodicity 
and thus the discreteness of both phase space variables , the area of the torus is characterized 
by a discrete Planck's constant, 

h = 2n/M (2.2) 

where M is the number of sites in both the coordinate and the momentum directions in 
phase space. Because of this, quantum dynamics defined by the cat map is considered to 
describe quantum resonance. Note that this is not a generic feature for other systems which 
have continuous phase space. 

The action S(Qj+i,Qj) for this linear map is easily constructed from conditions 

dS(Q j+1 ,Qj) . dS(Q j+1 ,Qj) 

dQ j+1 ~ Pj+h dQj ~ Pi (2 ' 3) 

Combining (2.1) and (2.3) gives 

S(Q j+1 , Qj) = ±(aQ 2 - 2Q j Q j+1 + dQ 2 j+1 ) (2.4) 

Before imposing the periodic boundary conditions, the propagator is 

U(Q j+1 , Q 3 ) = exp ^(aQ) - 2QjQ j+1 + dQ) +l ) (2.5) 

With periodic boundary conditions, one needs to sum over all equivalent initial points, thus 
yielding 

UiQf+uQj) = E ex P [^(a(Q J +27rm) 2 -2(Q j + 27rm)Q, +1 + dQ 2 3+1 )} 

= C(T,Ar)exp[—(aQ 2 J -2Q J Q J+l + dQ 2 +1 )] (2.6) 
where C(T,Af) is a constant depending on the form of T and M. 



7 



In fact , C (T, A/") vanishes in many choices of T and this sum gives a nontrivial value for 
the progagator only if the matrix has a special form. We choose 



for the elliptic case and 



r "= 3 2 ■ < 2 ' f 



for the hyperbolic case. 

For the special choice of the matrix elements Ti and T 2 made above, the propagator takes 
on the simple form, 

i 



UiU + U) = \l— r exp[—Q,Qi+i), (2-9) 



U2U + = y-^ex P [-(g| - Q,Q i+1 + Q 2 j+1 )]. (2.10) 

Since each iteration describes a permutation among sites, each site belongs to a periodic 
orbit. Thus quantum evolution follows the classical dynamics, resulting in the recurrence of 
the wave function (or equivalently, of the Wigner function) 

2.2 Decoherence in the Quantum Cat Map 

We now couple the system to a bath of iV harmonic oscillators linearly. []. The Hamiltonian 
of the bath of oscillators with coordinates q a and momentum p a (a = l,..,N) is 

N 2 2 2 

^ = E(y + ^)- (2-H) 

The interaction Hamiltonian between the system Q and the bath variables q a is assumed to 
be bilinear, 

N 

H C = Y, C -Qla- (2-12) 

a=l 

where C a is the coupling constant of the ath oscillator. 

By integrating out the bath variables, we get the reduced density matrix, 

PriQ^Q'pt) = J C =lC /g^ex P ^5(Q)+5 c (Q,g a )+^(g Q )-5(Q')-^(g',ga)-^(g;)]. 

(2.13) 

4 Nonlinearity in the coupling could enhance the interaction between the system and the bath. One mode 
in the system couples with numbers of different modes in the bath in the presence of nonlinear coupling. 
Then this type of coupling effectively increases the number of system in the bath N. Consequently, one can 
expect it helps decohere the system. Nonlinearity in the bath seems to have a similar effect. p5[ 
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where S is the classical action of the system defined in (2.4) and Sb, and Sc are the actions 
for the bath and the coupling respectively. The evolutionary operator J r for the reduced 
density matrix from time steps j to j + 1 is 

JriQj+hQ'j+i I Qj,Q'j,t) = J DQDQ' exp^-[S(Q) - S(Q') + A(Q,Q')] (2.14) 
in a path-integral representation [|26], [36], [H], [3£| , where 

l -A(Q,Q') = X> f ds [' dJr(s)[-ifi(s - s')R(s') - u(s - s')r(s')} (2.15) 
h h Jo Jo 

is the influence action. Here r = \{Q — Q'), R = \{Q + Q'), and /x(s), u(s) are the dissipation 
and noise kernels respectively ||38|| . 

If we consider the simplest case of an ohmic bath at high temperature kT > hA » hu a 
36| , and consider times shorter than the dissipation time, then we obtain a Gaussian form 



for the influence functional, with j^A(Q, Q') = — 2M J ! kT S 3 r|. where the noise kernel becomes 
local u(s) = 2M / ~fkT5(s) and 7 is the damping coefficient. Q The unit-time propagator 
becomes 

MQj+iiQj+i I QiiQ'j) = (Jr(Qj+i>Q , j +i I Q^Q'^O) = (exp^[5(Q J+1 ,Q J )-5(Q; +1 ,Q;)+2^ +1 ]). 

(2.16) 

Here ^ is a Gaussian white noise given by 

(0 = 0, (exp|er)=exp[-^^r 2 ] (2.17) 

where ( ) denotes statistical average over noise realization £. 
For the elliptic map, we get 

J r (Q j+ i,Q' j+ i I Q^Q'pO = (-^) 1/2 exp[j(- rj R j+1 - ■ 0-,.,) • (2.18) 

and for the hyperbolic map, 



J r (Q j+1 ,Q' j+l I Qj,Q'j,€) = {— r ) 1/2 exp[—(2r j R j + 2r j+1 R j+1 -r j R j+1 -r j+1 R j +£r j+1 )} 



(2.19) 

The Wigner function is defined as 

W(R,p) = ±J^4(R + r )p(R-r)e*p[=¥ L )dr. (2.20) 



Because the chaotic trajectory washes out information about the past rapidly, we expect that memory 
effect would be less important in classically chaotic systems than classically regular systems. Nevertheless, 
in the quantal case in which the classical stretching and folding behavior is suppressed, it would be still 
interesting to study how the non-Markovian behavior competes with nonlinearity pEf . We will discuss this 
issue in the next section. 
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where p is the momentum conjugate to r. The propagator Kt for the Wigner function is 



K T (R j+1 ,p j+1 | Rj,Pj,0 = ±rZ rj 12 rj+1 J r (Q j+1 ,Q' j+1 | Q j ,Q , j ,g)exp2L(p j r j -p j+1 r j+1 ). 

(2.21) 

This is reduced to the form of the classical cat map. For the elliptic case, 

R J = -p J+l + ^ Pj = Rj +1 . (2.22) 

For the hyperbolic case, 

Rj = 2R j+1 - p j+1 + f , pj = -3R j+1 + 2p j+1 - 2f . (2.23) 



Thus, without the environment, quantum evolution follows classical permutation |J3|. We 



can also say that the transformation from the classical map to the corresponding quantum 
propagator T — > Kx(Rj+i,Pj+i \ Rj,Pj,£ = 0) preserves the group structure. When coupled 
to the bath, the cat map is exposed to a Gaussian noise from the environment in each time 
step. The phase space is divided by a finite number of different periodic orbits and the period 
is known to increase roughly proportional to M . The discretized noise induces transition 
between different periodic orbits in an irregular way. As a consequence, the recurrence of 
some physical quantity will disappear in the quantum map and the classical type of mixing 
is regained. 

Fig.l shows Trp 2 , the linearized entropy (with the reversed sign) for various cases. If 
there is no interaction with the environment, the entropy is constant for both the regular 
and chaotic cases. When interaction sets in, Trp 2 r decays exponentially, showing that the 
system rapidly decoheres. There is no recurrence of this quantity observed. In spite of the 
discreteness of the points on the torus, we expect that the system behaves classically due to 
the influence of the environment. 

Note that the system also decoheres in a similar manner when the system is regular but 
with a slower rate |T7|] . These results indicate that if the underlying classical system shows 



chaotic behavior, even after quantized, the system still posseses the mixing behavior. This 
mixing property enhances the random perturbations from the environment, thus accelerating 
the suppression of quantum interference. However, in this particular example, the dynamics 
is essentially classical as is seen in (2.24) (In this case, the value of H comes in through the 
number of sites Af). More general cases should be examined. 

In Fig. 2, we show the mean displacement of points in the phase space from the initial 



configuration as a function of time. This is defined by I = J (Ax 2 + Ap 2 }, where Ax and Ap 
are the displacement from the initial phase space points, ( ) is the average over the phase 
space points and noise realizations. In the chaotic case, we see the recurrence disappears 
with just a small amount of noise (Fig.2a). Whereas in the regular case, the same amount 
of noise does not alter the qualitative picture of recurrence (Fig. 2b). In both cases, the 
decohered quantum system behaves close to the classical picture in which the regular and 
chaotic dynamics are clearly distinguished. For the elliptic map, the classical dynamics is 
completely periodic. For the choice of T\ in (2.7), the period is four. On the other hand, 
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for the hyperbolic map T2 in (2.8), the classical dynamics is nonperiodic. In spite of the 
discreteness of the points on the torus, the system behaves effectively classically due to the 
effect of the environment. 



3 Nonlinearity and Decoherence 
3.1 Quantum Kicked Rotor 

The kicked rotor and its map version, known as the standard map, are one of the most 
intensively studied models from both the quantum and classical points of view ||. The 
Hamiltonian of the kicked rotor is given by 

2 

H = J ti + KcosxE 7=-oot{t - J) (3-1) 

which describes a one-dimensional rotor subjected to a delta-functional periodic kick at 
t — j. Here x is the angle of the rotor with period 2tt, m is the momentum of inertia, p is an 
angular momentum, K is the strength of the kick and, in this case, the nonlinear parameter. 
When K > K c = 0.9716, the system becomes chaotic over the entire phase space. The 
average energy p 2 is known to show diffusive behavior like that of a Brownian particle under 
a stochastic force. This suggests the emergence of randomness in a deterministic chaotic 
system. 

The quantum dynamics of the kicked rotor is depicted by the corresponding Schrodinger 
equation 

d h 2 d 2 

ih—ip(x,t) = - 9 i>(x,t) + K cos xT,5(t - j)ip(x,t) (3.2) 
at 2m ox 

where ip is the wave function for the rotor. Denoting ipj as the wave function ip(x, t) at each 
discrete time t = j, and integrating (3.2) from j to j + 1, we obtain 

1>j+i(x) = ex PH^J^] exp[-z K( ^ SX ]^(x) (3.3) 

The quantum kicked rotor (QKR) is known to exhibit dynamical localization. After some 
relaxation time scale, the wave function becomes exponentially localized in the momentum 
space ||22|| . This may be interpreted as a particle moving in a lattice with a quasi-random 
potential. This heutistic picture seems to justify the analogy between the quantum kicked 
rotor to the tight binding model with an exponentially decaying hopping parameter which 
is known to show Anderson localization. The explicit transformation to the tight binding 
model was constructed in |23| . In spite of the nonrandom, deterministic nature of the kicked 



rotor Hamiltonian, numerical results there assuming sufficient quasi-randomness seem to 
support this analogy. Dynamical localization in this context arises from the suppression of 
classical diffusive behavior in the quantum dynamics. As shown by Ott et.al. [E7J], a small 



external noise can break the localization. Numerically they observed three different regimes 
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in the behavior of the diffusion constant D depending on the magnitude of noise. When the 
noise becomes sufficiently large, the system recovers the classical diffusive behavior. 

It is of interest to interpret the above results due to noise from the microscopic open 
system point of view. When the system interacts with an environment, we know that coarse 
graining of the environmental variables is also a source of noise and dissipation. We shall 
now derive the influence functional for the quantum kicked rotor and study its behavior. 



3.2 Quantum Kicked Rotor in a Bath 

Cohen and Fishman studied the case for the ohmic bath in detail ||25|| . They calculated 
explicitly the diffusion constant and the relevant time scales in terms of the noise correlation 
and the nonlinear parameter. 

We introduce a linear coupling of the system momentum p with each oscillator coordinate 
q a (a = 1, ..N) in the bath in the form H c = S^ =1 C a g Q p (Here q, p without the subscript 
a denote the system coordinate and momentum variables). For an ohmic environment, the 
action functional has the same form as (2.15), except that the coordinate variable Q is 
replaced by a momentum variable p. 

~A(p,p') = -5 ds ds'p~(s)[-ifj,(s - s')p + (s') - u(s - s')p~(s')} (3.4) 
n fi Jo Jo 

where p ± (s) = p(s) ±p'(s). 

In a similar way, we introduce the noise £(r) such that, 



r t s 

(exp[— % J £(t)p(t)]) = exp — — J ds J ds'p~(s)v(s — s')p~ 



(3.5) 



As before, we will examine processes in the time span where dissipation is small, thus ignoring 
the effect of the dissipation kernel /x(s). 

Under these assumptions, the action of the noise kernel can be formally absorbed in the 
propagator for the wave function. The unit time propagator for the wave function U^(J+1, j) 
is given by 

iK^ cos x zp^ ^^,P 

Ut(j + 1, j) = exp[ ] ex P[-7^;] ex P[ — H ( 3 - 6 ) 

where, as before, the noise term £ arises from using a Gaussian identity in the integral 
transform of the term involving the noise kernel in the influence functional |26|. Summing 
over all noise realizations ( ) gives the desired reduced density matrix, 

pAp, p') = (Vto(p) ^W)> (3-7) 

where 

^ +M (p) = U ( (j + l, ./) c Jt d>) (3.8) 

and if)j g(p) is the wave function under the influence of a particular noise history represented 
by£. ' 
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In the same way, Trp 2 T can be expressed as 



Trp 2 r = (EpE^tip) Pttf) MP') #(P))«' (3-9) 

where ( denote the statistical average of all possible noise histories of two independent 
noises £(t),£'(t) defined at each interval from j to j + 1. At high temperatures, £(t),£'(t) 
are reduced to two time-uncorrelated independent Gaussian white noises defined at each 
time step. 

There are many possible ways of introducing an interaction between the system and the 



bath, though many of them are related to each other as shown in [25]. One interesting case 
is when we intoroduce the coupling through the coodinate x. Then to preserve periodicity 
of the Hamiltonian under the coordinate transformation x — > x + 2ti, we need to restrict the 
range of noise to £ = nh(n = 0, ±1, ±2, ..), or choose the interaction Hamiltonian to be He = 
C a q a cos(x). In the latter case, we may further assume the form He = C a q a cos(x + (f) a ), 



where (f) a is the random phase [g7j to remove the coordinate dependence of the interaction. 
However, they all give the same result, but with different noise correlations. For example, 
from (3.6) we can calculate the propagator U^(j, 1) from t = 1 to t = j as 

1) = j - 1) U& -1,3- 2)...^(2, 1) 

= exp[-^f^] exp[-g] exp[-^] x .... (3.10) 

2 i 



'A coji l ovrl [_f!j CXp[ i Ml)P 

x exp[-^f^] exp[-g] exp[-^£ 
exp[- iKco8 ^ j)) ] exp[-g] x ... x exp[- ^ cos( ^ (1)) ] exp[-g^ 



where = £(j) + ... + £(1). Thus, this describes the same dynamics as couplings through 
x via He = C a q a sinx, as long as the noise r](j) remains small. The correlation of the two 
different noises rj(j) and are related by 

( V (rHr')) = EUEiumm) (3.H) 

Note that even if £ is a white noise, r\ is not necessarily white. 

3.3 Localization and Decoherence in the Quantum Kicked Rotor 

The eigenvalue equation for the quasi-energy state in the QKR is given by [p3[ 



2 

exp(--K cosx) exp(--—)u UJ (x) = exp(- -w)u u (x) . (3.12) 

This can be transformed to 

{r i / V 2 mt .Kcosx^ r .Kcosx.A . . i ,.1 

6XP[ "^ ( 2^ ~ ^ )][ ~ 1 tan( ^ _ )] ~ [ + 1 tan( ^~ )] j [ + ex P("^ K cosx )] 2 = °' 

(3.13) 
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If we define lit,; [ X ) cLS 



i 1 

u u (x) = [1 + exp(— -K cosx)]-Uu(x) (3.14) 



then (3.15) can be written as 



[l + expK(£-u,))]-tan(^) 
Expanding ^(x) as w w (x) = T lk z™ 00 Uke' lkx , we get 



m w (x) = 0. (3.15) 



where 



T k u k + T, r ^ W r Uk+ r = Eu k (3.16) 



T k = tan(— — ) (3.17) 

W k = -L p Me™tanfi^) (3.18) 



and E = —Wo- As pointed out in p3| , (3.16) gives the eigenvalue equation for the tight bind- 
ing model which describes electron motion in a quasi-random potential T k with the hopping 
parameter W k . The property of this model depends on the rationality of the coefficient of k 
in the potential. When the coefficient is irrational, the model is known to show Anderson 
localization. 

We can now analyze the relation between the breaking of dynamical localization and 
quantum decoherence. Loss of quantum coherence is measured by the density matrix be- 
coming approximately diagonal. Decoherence in the quantum Brownian model has been 
studied extensively for this problem. We refer the reader to recent work on this topic ||16|| . 
In Fig. 3a we plot the linearized entropy Trp% versus the energy (p 2 /2) in each diagram (we 
set the mass m = 1 for all the numerical calculations). Note that the two effects are cor- 
related to each other as expected. This shows that derealization occurs as the quantum 
coherence breaks down, suggesting that derealization and decoherence occur by the same 
mechanism. As the noise strength increases, we see that decoherence works more efficiently. 
Also as the nonlinearity parameter K, the strength of the kick in this model, increases, the 
system decoheres more rapidly (Fig.3b). At the same time, the amount of derealization 
measured by the diffusion constant also increases. For all the numerical results presented 
in this paper, we use the Gaussian wave packet as the initial condition. However, we also 
checked that the qualitative results given in this paper are insensitive to the initial condi- 
tion. This may be viewed as another characteristic of the chaotic system. For QKR, as the 
nonlinear parameter K decreases, the results become more sensitive to the initial condition 
28fl . In nonchaotic systems, the sensitivity in this sense can be used to choose a preferred 



initial state in accordance to some specific criterion, such as least entropy production [p9 
etc. (For chaotic cases, see ]30 for a related argument but from a different point of view). 



This may be explained in the following way. Because we use a coupling through the 



momentum, the time scale for the system to lose coherence is given by td e 



7 V 8p 
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where Xub — h / \/2irmkT is the thermal de Broglie wavelength, and dp is the relevant 
momentum scale. After this time, noise will destroy the quantum coherence between these 
momentum separation. In the kicked rotor case, localization will occur due to the coherence 
around 5p ~ A, where A ~ lh is the localization length. Since I ~ K 2 , this gives t^ e 
Therefore, nonlinearity increases the rate of decoherence. 



The relation between the diffusion constant D and the noise strength is given in |27|, [25 
For our case, K/h ^> 1 and for weak noises, we can consider the particle as undergoing a 
random walk with hopping parameter l/t c . Then D = A—, where A gives the localization 

''dec 

length. From this, we get 

£> = 7(4-) (3-19) 

A 0dB 

4 Decoherence and Irreversibility 

The Wigner function representation is often used to examine the quantum to classical tran- 
sition. For a linear system the Wigner function is known to show a smooth convergence 
to the classical Liouville distribution. But if the system Hamiltonian has a nonlinear term, 
quantum corrections associated with the higher derivatives of the potential pick up the rapid 
oscillations in the Wigner function and it no longer has a smooth classical limit | 40|| . How- 



ever, upon interaction with an environment, a coarse-grained Wigner function can have a 



smooth classical limit [f|l| for nonlinear systems 



The Wigner function at time t — j is defined as 

1 +2n 11 
W s {X,p) = — J dy Pj (X ~^y,X + -y). (4.1) 

-2tt 

where X = + x'), y = x — x' . From (3.3), the unit time propagator for the Wigner 
function of QKR is, 

W j+1 (X,p) = exp[-^^A p ] expi-^-dxjW^X^) (4.2) 

where A p = exp[|<9 p ] — exp[— ^d p ) measures the effect of the kick. We can see the effect of 
quantum corrections more clearly if we expand A p in orders of h: 

A p = exp[p p ] - exp[-f«9 p ] = %d p + ^ + ^- Q d 5 p + (4.3) 

With this, the first exponential in (4.2) contains the classical propagator times quantum 
corrections of even orders of K. 

exp[-^fAA p ] = exp[-K smXd p ]exp[-^K sinXd^]exp[-^K sinXd 5 p ] (4.4) 

If the initial system wavefunction is discribed by a Gaussian wave packet with width 
Sp(» h), , we would expect to see a classical-like evolution of the packet at short times. 
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When the width of the contracting wave packet gets small, and becomes comparable to 
H, the effect of quantum corrections appears, namely, the corrections from the exponent 
which is higher order in h in (4.4). By comparing the classical and the quantum terms, 
we can easily evaluate the length scale at which quantum corrections become important,i.e. 
when 5p(t) ~ h. Here Sp(t) = 5p(0)e~ xt , where A is the Lyapunov exponent given by 
A ~ ln(K/2) . As shown in (1.2), from this expression, we can define the time scale tg as 
t E ~ ln Sj ^/X. Because the Wi gner function or the expectation value of any observable 
follows classical trajectories when t < t E , this has been called the Ehrenfest time. Note that 
in the continuum case, this definition gives us a different time scale for each term in the 
expansion [Q. Hereafter, we set m = 1 for brevity. 

If the interaction with the environment has a form in (3.6), the major effect of the bath 
is the appearance of a diffusion term in (4.2), such that, 

W j+l (X,p) = exp[D x d 2 x ] exp[- Ks ^ X A p ] ex V [-pd x }W 3 {X )P ) 

ft 3 

« exp[D x d 2 x ] exp[~—KfanX6%\W j (X - p + KsmX,p - K sin X) (4.5) 

where D x = IM^kTTi is related to a constant of the noise kernel v{s) defined before (2.16). 

Competition amongst the three terms with different physical origins is apparent: The 
first term in (4.5) is the quantum diffusion term, the second is the quantum correction term, 



and the third is purely classical evolution. As discussed by Zurek and Paz if D is 
sufficiently large, the effect of quantum corrections becomes inconspicuous. In this case, the 
diffusion term traces out a small scale oscillating behavior before quantum corrections have 
a chance to change the classical evolution. Then one may expect the time evolution of the 
Wigner function to be like that of classical evolution with noise. In this case, we can ignore 
the quantum correction in (4.5) and write the evolution equation as 

W j+1 (X,p) = exp[D x d^]Wj(X-p + KmnX,p-KwnX) (4.6) 

The role of quantum diffusion is to add some Gaussian averaging so that the contracting 
direction in phase space will be suppressed while it does not affect the stretching direction. 
As long as the width of the wave packet is large such that the first term is negligible, the 
evolution should be Liouvillian (time reversible if we assume infinite measurement precision). 
Furthermore, we expect that after the width of the packet along the contracting direction 
becomes comparable to the diffusion generated width (in the Gaussian wave packet), the 
dynamics will start showing irreversible behavior arising from coarse graining (as distict 
from irrreversibility from instability). Consequently, entropy should increase in this regime. 
In Fig. 4, we plot the von Neumann entropy for the dynamics of (4.6). We can see three 
qualitatively different regimes: I. the Liouville regime: the entropy is constant and the dy- 
namics is time reversible. II. the decohering regime: the entropy keeps increasing due to 
coarse graining. III. the finite size regime: due to the bounded nature of the phase space, 
the entropy shows saturation. Our result from quantitative analysis seems to confirm the 
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qualitative description of Zurek and Paz [[18] who used the inverted harmonic oscillator po- 
tential as a generic source of instability. Since the phase space in their model is not bounded 
they do not see Regime III. Similar features appear in the quantum cap map (Fig. 5). In this 
case, the full quantum dynamics can be calculated in a simple way. Resemblance with the 
result of a classical rotor with noise is obvious. However, in this case, the stable entropy is 
smaller than the maximum value which may be explained as a finite (phase space) size effect. 



Quantum diffusion defined by the spreading of the wavefuction is known to be dynam- 
ically stable. The authors of |3l| performed the time reversal at some time and saw the 
diffusion constant and even the wave function itself coming back to the same state within 
the accuracy of computation. As we know, these time reversal behaviors cannot be seen in 
the classical case due to the instability of the trajectories. This is also true in real physical 
systems for which one can access information with only finite precision 0]. 

In Fig. 6, we perform the time reversal at t = 200. In the quantum case without bath, 
the system completely returns to the original state after exactly the same amount of time. 
Thus, the system is highly stable in spite of its random appearance. On the other hand, 
in the classical case, instability prevents the reversibility even without interaction with the 
bath. When the interaction is turned on, we see that the reversibility in the quantum system 
is gradually lost, and irreversibility appears as the noise strength increases. 

In a realistic physical system which has a finite precision due to numerical or instrumen- 
tal limitation, we expect this type of irreversibity is inevitable for the chaotic system even 
without noise. If the minimum precision in length is denoted as e, the time scale up to which 
the deterministic picture is valid is determined as t p ~ -log— . At t > t p , the system starts 
losing information about the past. This may be the source of irreversibility for the classical 
chaotic system. For the system we are studying, t p » t E > 1 holds. Then the quantum 
effect smears the contracting evolution of the region in the phase space before information 
about the past is lost. Therefore, the wave packet traces back the same trajectory as it 
evolved. If t c < t p , we see the irreversibility from coarse graining before the limitation of 
the measurement becomes evident. This type of irreversibility is another characteristic of 
classicality peculiar to chaotic systems. 
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Figure 1 The linearized entropy (with opposite sign) is plotted against time. If there is no 
environment, the entropy is constant for both the hyperbolic and elliptic cases, indicating the 
purity of the state (dotted line). For the hyperbolic map, even though classically this system 
is strongly chaotic, the corresponding quantum system does not show chaotic behavior. This 
situation changes drastically when the system interacts with a thermal bath. In this case, 
the entropy due to coarse graining keeps increasing. Note that in the hyperbolic case (solid 
line), the rate of entropy increase is larger than in the elliptic case (dashed line). N = 50 is 
used here (also in Fig. 2). 

Figure 2 The mean phase space point displacement is shown. When there is no environment 
(dotted line), the system shows recurrence in both the hyperbolic (a) and elliptic (b) cases. In 
the presence of an environment (solid line), the hyperbolic map loses the recurrence behavior 
under the Gaussian noise with a = 0.08 and maintains a near-constant value, indicating 
ergodicity of the classical hyperbolic map. On the other hand, the elliptic map still shows 
recurrence with the same amount of noise, corresponding to the classical periodicity. 
Figure 3 Trp 2 (solid line, left scale) and < p 2 /2 > (dashed line, right scale) are plotted in 
Fig. 3a versus time for K = 10 and H — 1. Three noise strength values a = 0, a = 1.0, 
a = 2.0 are plotted here, corresponding to the family of lines from up to down for Trp 2 
and down to up for < p 2 /2 > (note that Trp 2 = 1 for a = 0). As the noise strength 
increases, the decoherence time shortens, and Trp 2 . decays rapidly. This accompanies the 
increase of diffusive behavior in < p 2 /2 >. In Fig. 3b, the same observables are plotted but 
with different i^-values. The upper, middle, and lower solid lines (lower, middle, and upper 
dashed lines) correspond to K = 0.5, K = 5, K = 10, respectively. Here a = 1.0, % — 1 are 
fixed. We see that increasing nonlinearity shortens the decoherence time. Note that when 
K = 0.5, the system merely diffuses. 

Figure 4 The von-Neumann entropy is plotted versus time for the kicked rotor. For the first 
25 steps, the system does not produce any entropy. The evolution is reversible. Transition 
sets in in the next few steps, the dynamics changes its character from reversible to irreversible. 
Because the nonlinear term is from the sinusoidal function, the onset of this regime is different 
at every phase space point. Consequently we can only see the averaged behavior through the 
entropy function. Around t — 40, saturation occurs due to the finiteness and the periodic 
nature of the phase space (a = 0.1 for this case). 

Figure 5 The von-Neumann entropy is plotted versus time for the quantum cat map. Due 
to the simplicity of the system, we see the same qualitative features as in Fig. 4. The 
entropy starts increasing around t = 55 and maintains a near-constant rate of production 
while showing large oscillations. After 200 steps, entropy production seems to saturate and 
starts decreasing into some stable value. Af = 64, a = 0.04 are used here. 
Figure 6 The time reversal is performed at time = 200. For QKR without a bath (lower 
curve), the system posseses time reversal invariance f3"I| . When the interaction with the 
environment increases, the system gradually regains irreversibility as observed in classical 
systems. Here the noise strength a = 0.5 (middle curve) and a = 1.0 (upper curve). Also 
h = l, K = 10. 
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